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Introduction 

Strongly Interacting Massive Particle (SIMPS), by which we will always mean 
neutral, stable SIMPs, are of current interest for at least three reasons: 

• They could be a dark matter constituent as suggested some time ago by Dover, 
Gaisser and Steigman [1] and by Wolfram [2]. Starkman et al.,[3] show SIMPs 
would be restricted to rather narrow mass ranges if they were to exhaust Q — 1. 
We will not make this assumption and will consider SIMPs outside the regions 
allowed by the analysis of ref . [3] . 

• It is possible that the lightest SUSY particle (LSP) is strongly interacting and 
hence, if R-parity is conserved, would form a colorless SIMP. Possibilities, such 
as a gg bound state are discussed in ref. [4]. 

• An explanation of the ultra high energy cosmic ray events (UHECRs) proposed 
by Farrar, Kolb and co-workers [5] is that they are due to interactions of SIMPs 



* Presented by Vigdor L. Teplitz. 

^To appear in the proceedings of the International Conference On Orbis Scientiae 1999: Quantum 
Gravity, Generalized Theory Of Gravitation And Superstring Theory Based Unification, 16-19 De- 
cember 1999, Fort Launderdale, Florida. 
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with a mass below 50 GeV and a cross section for interactions with nucleons on 
the order of a (few) millibarns. 

This summary will review two laboratory experiments that might detect SIMPs. For 
more detail, the reader is encouraged to examine with care ref. [6] and the paper on 
which it is based [7]. 

In Section 2 we consider the possibility of finding SIMPs bound in ordinary nuclei 
by searching for anomalously heavy isotopes of high-Z nuclei. It is a pleasure to note 
that the accelerator mass spectrometry (AMS) group at Purdue is in the process of 
performing the experiment^ suggested in ref. [6]. In Section 3, we address the extent 
to which production and detection of SIMP-anti-SIMP (SS) pairs might be performed 
at the Tevatron. 

Our results, in brief, are that the AMS experiment should be sensitive to SIMPs 
over a wide range of parameter space: (ctsn,Ms), where Ms is the SIMP mass and 
<jsn is its cross section for scattering off nucleons. The Tevatron, on the other hand 
is likely only to produce and to detect SIMPs in a much more restricted range, but 
one that includes much of the mass range for which the SIMP could be the UHECR 
explanation. It would be only fitting, since much of the work on that possibility [5] was 
done at Fermilab, if SIMPs were to be detected at Fermilab and we encourage those 
with influence in the collaborations to explore vigorously that possibility. Finally, we 
note that we proceed without committing to a specific SIMP model. We parameterize 
the experimental predictions in terms of the two parameters ctsn and Ms- 

SIMPs in Nuclei 

We know a fair amount about SIMP binding in nuclei from the phenomenology 
of hyper-fragments. See, for example, Povh [8] for a readable review. Based on that 
experience, we can write for the binding B of the SIMP in a nucleus A the relation: 

B = \V sn \-t: 2 /(2hR 2 ) , (1) 

where \i is the reduced mass of the S-A system, R is the radius of the nucleus A, and 
V is the S-N potential averaged over the volume of the nucleus X. We expect the low 
energy potential, Vsn, to be always attractive. This is true if exchange of vacuum 
quantum numbers dominates. We assume this to be the case, and have not found a 
model to the contrary. Under this assumption, the SIMP can be bound in a nucleus 
for which ft and R 2 ~ A 2 / 3 are large enough to make the kinetic energy less than the 
(average) magnitude of the attractive potential. 

From equation (1) we see that the best chance of finding SIMPs is to search 
in high Z (large) nuclei which minimize the kinetic energy term. Capture by light 

$We are grateful to Professor Ephraim Fischbach for keeping us informed as to progress on this crucial 
experiment. 
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elements at the time of cosmic nucleosynthesis has been studied in ref. [9]. Atomic Mass 
Spectrometer (AMS) searches to date are reviewed in the careful study of Hemmick 
et al., [10] where one learns the somewhat surprising fact that previous searches have 
only been conducted up to sodium (Z = 11, A = 23). This makes the current Purdue 
AMS experiment particularly exciting. They are looking in gold (Z = 79, A = 100). 
At the risk, however, of sounding greedy, we would be interested in anyone with AMS 
equipment, a supply of Lawrencium (Z = 103, A = 262), and an affinity for maximal 
experiments. However, as will become clear below, we do not want any such calls from 
interested parties to be collect. 

How big is the potential Vxjv? We take this as a parameter, but we can put 
an approximate LOWER bound on it from the requirement that primordial S and S, 
left over from the early universe, not overdose the universe so that it couldn't have 
continued expanding until today (late 1999). The classic book of Kolb and Turner [11] 
tells us that the number density of primordial SIMPs behaves as 

ns-iMsass)- 1 . (2) 

Equation (2) says that too small an annihilation cross section means too many 
SIMPs will be left over from the early universe, and Kolb and Turner collect together 
the numerical recipes for computing how small is too small. We still need, however, to 
relate the annihilation cross section, a ss - to the SIMP-nucleon cross section, <tsn and 
to the S — N potential in Equation (1). We make the simple ansatz 

Vsn = Vnn {&sn/(Tnn) 1/2 (3) 

°SN = f3°NN V S S ( 4 ) 

where (3 should be on the order of one. Note that V S n goes as f3 1 ^ so that our results 
for binding will not be highly dependent on the precision of Equation (3). 

Now that we know, for each point in the Ms, <Jsn parameter space, the primordial 
S abundance and the binding energy in nuclei, we are almost ready to compute for our 
friends at Purdue, the abundance of anomalous gold-gold with a SIMP bound in it. 
First, however, we need a scenario for how the SIMPs get bound into the gold. Our 
picture is as follows: 

• We assume that the ratio of SIMPs to protons in the galaxy is the same as the 
cosmic ratio, but that most of the SIMPs are in the galactic halo (i.e., that their 
density distribution is p ~ R 2 ), not in stars. We can then calculate the SIMP 
flux on the Earth, since we know that the Earth is traveling through the galaxy 
with a velocity of about 200km/ s which not too different from the galactic virial 
velocity. 

• We assume that when the SIMP hits the Earth, it is slowed by scattering with all 
nucleons and nuclei at a rate determined by as n, but can only be captured by a 
nucleus that is large enough. 



3 



• Gold must compete, for SIMP capture, with the most abundant nuclei large 
enough to bind the SIMP. Our comparative estimates use, as the most abundant 
elements:, aluminum (A = 27), barium (A = 137), and lead (A = 206). 

Our procedure is then as follows: 

• We chose values for Ms and <tsn and then determine whether, for that point in 
parameter space, there is binding in gold. 

• Assuming that there is binding, we then determine (a) the mean free path in 
Earth from the galactic virial velocity and ctsn, and (b) which of the 3 elements 
above is gold's chief competitor for SIMP capture. 

• From the ratio of the abundance of gold to its chief competitor, the mean free 
path, and the average density of Earth, we then compute the chance of a particular 
gold nucleus within a mean free path to capture an incident SIMP. Multiplying 
by the flux (see above) of SIMPs and the time for which the sample being put in 
the AMS target has been exposed (which is why we don't want collect calls from 
those long in Lawrencium) gives us the fraction of gold nuclei in the sample that 
should have a SIMP if they exist at that point in parameter space. 

Finally, we assume^ that the exposure time is 10 million years because there are re- 
gions that are geologically inactive over such periods and have had for example "placer" 
gold in the beds of streams for a longer period than thatrl 

The results are shown in the table. It gives log 10 of the ratio of normal to anomalous 
gold nuclei. The dashes indicate parameter values for which there is either no binding 
in gold or overclosure of the universe. One sees that smaller values of <tsn give larger 
ratios of anomalous to normal gold. This is because smaller values imply that only lead 
has a nucleus large enough to compete with gold for SIMP capture and because the 
smaller cross section means more primordial abundance. The important thing to take 
away from hours of table study is the fact that the relative abundance entries are all 
considerably higher (for anomalous to normal) than the limits of 10~ 20 that have been 
set in AMS work on some of the light elements. This provides reason to expect that, if 
the SIMPs are there, the Purdue AMS people will find them. 

SIMPs at Fermilab 

Next we consider SS production at the Tevatron. Since we are talking neutral 
SIMPs, we expect little or no signal in the central tracker and in the electromagnetic 
calorimeter. However, in the hadron calorimeter, we expect to detect SIMP signals if 

§We appreciate conversations with Professor E. T. Herrin on searching for old exposed gold. 
^Purdue is, we believe, taking collect calls from people who have pieces of gold with such provenance; 
they give all but a small fraction of a mole back at the end of the run. 
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Table 1. M x (vertical) is in units of GeV, and a X N (horizontal) is in units of mb. 
Table entries are log 10 (l//), and the — indicates those cases for which X does not 
bind at all. 



(Jsn is large enough. The detection of SIMPs is possible if one triggers on two back-to- 
back hadron calorimeter showers, accompanied by little else. We will use 10 GeV for 
the minimum size showers for which such triggering might be done. Our task now is 
to determine: 

• For what values of {M S) <Js} will the SIMP interact in the steel plates of the 
hadron calorimeter? 

• For what values of these parameters will we get calorimeter showers greater than 
10 GeV or more? 

• Can one recognize a SIMP shower if one sees one? 

• How many such events should we expect? 

First we look at the region of parameter space for which there will be interaction. 
The minimum annihilation cross section permitted by the cosmology argument is 
~ 3 x 10~ 13 barns, which corresponds through Equation (3) to about a microbarn for 
the S-N cross section. SIMPs with such small cross sections won't shower in 1 meter 
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Figure 1. Contours for energy loss as a function of {Mx,Ex}- The contours displayed are in steps 
of 10 GeV. 

of steel, but for a higher cross section of a few millibarns, we would expect 10 or more 
interactions with the 10 27 nucleons/cm 2 in the 1 meter. 

To estimate the energy we expect in a shower resulting from a SIMP interaction 
in the steel plates of a hadron calorimeter we use a cosmic ray rule of thumb kindly 
provided by G. Yodhfl] who says that, in a high energy strong interaction, about half 
the center of mass energy goes into inelasticity. In the figure, we give the (laboratory) 
energy released in the calorimeter as we vary the mass and energy of the SIMP; the 
straight lines are constant shower energies. One sees that the bigger the SIMP lab 
energy, the greater a SIMP mass will result in a given shower energy. 

Consider now the question of whether we would recognize a SIMP shower if we saw 
one. The background for SIMP showers would likely be neutron showers and K decays. 
The distinguishing feature would be shower opening angle. A pion moving transverse 
in the cm. system would have a lab angle given by tan# = 1/7. Comparing the angle 
for a SIMP with that from a neutron of the same energy, the SIMP shower should be 
wider by roughly the ratio of the masses. 

Finally, we turn to the number of SIMP pairs the Tevatron might produce. We 
scale the (known) production rate of jets by the ratio of the S-N cross section to that 
of Meson-N, which we take to be on the order of 30 millibarns. So long as the SIMP 
energy is a few times its mass, we don't worry about phase space suppression. We 
assume conservatively a cross section of about 3pb for any one parton in the region 

"This useful approximation from Professor Gaurang Yodh made the whole trip to Paris (where the 
conversation took place) well worthwhile (and the food was OK too). 
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E > 200GeV. This implies about 6000 events in Run II. The estimate of [5] is that 
the Nucleon-UHECR cross section needs to be over a tenth the Meson-Nucleon cross 
section, so we estimate 600 events in the Tevatron run if SIMPs are the explanation for 
the UHECR events. 



Summary 

For the Table we see that there is SIMP binding in gold for MgUsN > 5mbGeV 2 , 
and that AMS experiments sensitive to one part in 10 20 can detect the existence of 
SIMPs of mass less than a TeV, while the region of interest for explaining UHECRs 
can be explored with a sensitivity of one part in 10 16 or less. Looking for SIMPs at the 
Tevatron is more difficult, but over half the region of interest for explaining UHECRs 
could be searched in the upcoming Run II by looking for (wide) back to back jets with 
no signal in the central tracker or EM calorimeter. Interested people are urged to send 
any gold in their possession for which they can prove 10 7 (or more) years of exposure 
to Purdue. Almost all nuclei will be returned intactQ. 

We thank D. Berley, K. Brockett, K. De, D. Dicus, M.A. Doncheski, R. Ellsworth, 
G. Farrar, E. T. Herrin, D. Rosenbaum, R. Scalise, and G. Yodh. The work of RNM 
has been supported by the National Science Foundation grant under no. PHY-9802551. 
The work of Olness, Stroynowski, and Teplitz is supported by DOE. 
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** Sorry, but neither Purdue, SMU, or Maryland can be responsible for nuclei lost or damaged in transit 
(by either AMS or Fedex) 
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